Influence of the contacts on the conductance of interacting quantum wires 
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We investigate how the conductance G through a clean interacting quantum wire is aflected by 
the presence of contacts and noninteracting leads. The contacts are defined by a vanishing two- 
particle interaction to the left and a finite repulsive interaction to the right or vice versa. No 
additional single-particle scattering terms (impurities) are added. We first use bosonization and the 
local Luttinger liquid picture and show that within this approach G is determined by the properties 
of the leads regardless of the details of the spatial variation of the Luttinger liquid parameters. 
This generalizes earlier results obtained for step-like variations. In particular, no single-particle 
backscattering is generated at the contacts. We then study a microscopic model applying the 
functional renormalization group and show that the spatial variation of the interaction produces 
single-particle backscattering, which in turn leads to a reduced conductance. We investigate how 
the smoothness of the contacts affects G and show that for decreasing energy scale its deviation 
from the unitary limit follows a power law with the same exponent as obtained for a system with a 
weak single-particle impurity placed in the contact region of the interacting wire and the leads. 

PACS numbers: 71.10.Pm, 72.10.-d, 73.21.Hb 



I. INTRODUCTION 

The low-energy physics of one-dimensional (Id) metals 
is not described by the Fermi liquid theory if two-particle 
interactions are taken into account. Such systems fall 
into the Luttinger liquid (LL) universality class^ that is 
characterized by power-law scaling of a variety of correla- 
tion functions and a vanishing quasi-particle weight. For 
spin-rotational invariant interactions and spinlcss mod- 
els, on which we focus here, the exponents of the differ- 
ent correlation functions can be parametrized by a single 
number, the interaction dependent LL parameter K < 1 
(for repulsive interactions; K ^ \ in the noninteract- 
ing case). As a second independent LL parameter one 
can take the velocity Vc of charge excitations (for details 
see below). Instead of being quasi-particles the low ly- 
ing excitations of LLs are collective density excitations. 
This implies that impurities or more generally inhomo- 
geneities have a dramatic effect on the physical properties 
of LLs,2iii^ 

In the presence of only a single impurity on asymptot- 
ically small energy scales observables behave as if the Id 
system was cut in two halfs at the position of the im- 
purity, with open boundary conditions at the end points 
(open chain fixed point) ^i^ In particular, for a weak 
impurity and decreasing energy scale s the deviation of 
the linear conductance G from the impurity-free value 
scales as (s/so)^'^^^', with K being the scaling dimen- 
sion of the perfect chain fixed point and sq a charac- 
teristic energy scale (e.g. the band width). This holds 

as long as |T4ack/so|^ (s/sq)^'^"^^ < 1, with Vback be- 
ing a measure for the strength of the 2kF backscatter- 
ing of the impurity and kp the Fermi momentum. For 
smaller energy scales or larger bare impurity backscatter- 
ing this behavior crosses over to another power-law scal- 
ing G{s) ^ (s/sq)^*-^/^"^', with the scaling dimension of 
the open chain fixed point l/K. This scenario was veri- 



fied for infinite LLsSi^ as well as finite LLs connected to 
Fermi liquid leads r*^ a setup that is closer to systems 
that can be realized in experiments. In the latter case 
the scaling holds as long as the contacts are modeled to 
be "perfect" , that is free of any bare and effective single- 
particle backscattering, and the impurity is placed in the 
bulk of the interacting quantum wire. For an impurity 
placed close to perfect contacts the exponents change to 
2{K — 1)/{K + 1) (close to the perfect chain fixed point) 
and l/K — 1 (close to the open chain fixed point)i2ii£ 

The role of an inhomogeneous two-particle interaction, 
that is an interaction that depends not only on the rel- 
ative distance of the two particles, but also the center 
of mass is less well understood. In the present publi- 
cation we will fill this gap. Such an inhomogeneity will 
generically appear close to the interface of the interacting 
quantum wire and the leads and a detailed understand- 
ing is thus essential for the interpretation of transport ex- 
periments on quasi-ld quantum wires^ We here use two 
models to study the effect of two-particle inhomogeneities 
on the linear conductance. We first investigate the so- 
called local Luttinger liquid (LLL) modelfi2iiii4 that is 
characterized by a spatial dependence of the LL parame- 
ters K and Vc, with K = I and Vc = vp in the leads {vp is 
the noninteracting Fermi velocity) . We show that regard- 
less of the details of the spatial variation the conductance 
always takes the perfect value l/(27r) (in units such that 
h = 1 and the electron charge e = 1). Thus the LLL 
description cannot produce any effective single-particle 
backscattering from the contact region generated by an 
inhomogeneous two-particle interaction. Our results gen- 
eralize earlier findings obtained for step- like variations of 
the LL parametersii2iiii4 

We then study a microscopic lattice model with a spa- 
tially dependent nearest-neighbor interaction. Across the 
contact between the left lead and the wire the interac- 
tion is turned on from zero to a bulk value U and corre- 



spondingly turned off close to the right contact. We show 
that this two-particle inhomogeneity generically leads to 
an effective single-particle backscattering and a reduced 
conductance. To compute the latter we use an approx- 
imation scheme that is based on the functional renor- 
malization group (fRG)pifi We numerically and analyti- 
cally investigate the dependence of the conductance on U, 
the length of the interacting wire N, and the "smooth- 
ness" with which the interaction is turned on and off. 
For weak effective inhomogeneities we analytically show 
that l/(27r) — G displays scaling with the energy scale 
Sn = ttvf/N set by the length of the wire. The expo- 
nent we find is consistent with the one found for a weak 
single-particle impurity close to a perfect contact i2iiii We 
give the energy scale up to which this scaling holds. It 
depends on the bulk interaction U and the 2fci? Fourier 
component of the function with which the interaction is 
varied close to the contacts. The latter provides a quanti- 
tative measure of the "smoothness" of the contacts. This 
finding suggests a similarity between the universal low- 
energy properties of a LL with a single-particle inhomo- 
geneity and a two-particle inhomogeneity. We discuss 
this similarity but also point out differences. 

This paper is organized as follows. In Scc.|n]we study 
the two-particle inhomogeneity within the LLL picture. 
To motivate the LLL model in Sec. Ill Al we present a 
brief introduction into bosonization and the Tomonaga- 
Luttinger model. The transport properties of the LLL 
model are investigated in Sec. Ill Bl In Sec. lIIII we intro- 
duce our microscopic model and give the basic equations 
of the fRG approach. We then discuss our numerical re- 
sults in Sec, nil Al and the analytical findings in Sec. lIIIBl 
We close in Sec. II VI with a summary. 



II. THE LOCAL LUTTINGER LIQUID 
DESCRIPTION 



A. Model and generalized wave equations 

In order to also elucidate the limitations of the LLL 
picture we first shortly recall the basic ideas necessary to 
justify the Tomonaga-Luttinger-modeU*iSiii for interact- 
ing fermions in one dimension in the homogeneous case. 
We consider a system of length L with periodic boundary 
conditions. In second quantization the two-body interac- 
tion reads 



V = - 



v{x — x')5 p{x)6 p(x' )dxdx' 



+^AA2^(0) - -v{Q)N, 



(1) 



where 5p{x) = il)'^ {x)ip{x) ~ M jL is the operator of the 
particle density relative to its homogeneous value, with 
M the particle number operator. The Fourier transform 
of the two-body interaction is denoted as {J(fc). The last 
term is usually dropped as it only modifies the chem- 
ical potential. If the range of the two-body interac- 
tion is much larger than the mean particle distance only 
particle-hole pairs in the vicinity of the two Fermi points 
are present in the ground state and the eigenstates with 
low excitation energy^ This allows to linearize the dis- 
persion around the two Fermi points and to introduce two 
independent types of fermions, the right- and left-movers 
with particle density operators 8p± (x) — p± (x) — A/± 



r / \ \ ^ ik X 

dpa(x) = y2^e " p„ „ = — 
dx 



27r ^^ 

n=£0 



AkjiX 



~ Pn,a 



(2) 



As shown below the field operators <&« are convenient ob- 
jects for the solution of the problem. ^^ In the subspacc of 
low energy states the Fourier components of the density 
Pn.a obey the commutation relationsi^ 



[Pni.aiPn.p] = amSa/S S„ 



(3) 



In this section we discuss how sufficiently smooth con- 
tacts (perfect contacts) can be properly described in 
a LLL picture. This approach was successfully used 
to show that for perfect contacts the properties of the 
leads and not the finite size quantum wire determine 
the conduct anceii2iiii4 Up to now the LL parameters 
were always assumed to vary step-like. Here we present 
a simple derivation of this result for an arbitrary varia- 
tion of the LL parameters. Also the role of impurities 
in the interacting wire was investigated within the LLL 
pictureAi^ It was shown that the exponent of the tem- 
perature dependence of the conductance is different for 
impurities placed in the bulk and close to the contactsiS. 
Using the fRG this was later confirmed in a microscopic 
model in which the contacts were modeled to be arbitrar- 
ily smooth»i£. 



After proper normalization they take the form of boson 
commutation relations. 

One can write the operator of the kinetic energy as a 
quadratic form of the p± i^ The Tomonaga model then 
results from replacing 6p{x) by Sp+{x) + Sp-{x) and Af 
by Af+ +N- in Eq. iQJ. With the well known "g-ology" 
generalizationi the boson part of the interaction reads 



V. 



TL,f) 



{2g2{x - x')5p+{x)6p^{x') 



(4) 



+gi{x — x') [5p^{x)5p^{x') -\- 5p^{x)5p-{x')W dxdx' . 

This reduces to Tomonaga's original model for §2 = 
g4 = w, while Luttinger^ ■ later independently studied 
the model with 54 = 0. He also discussed the special 
case g^ix) ~ 5{x)^ which corresponds to the interaction 
term in the massless Thirring model. ^^ This simplifies 



various aspects but brings in infinities which have to be 
removed e.g. by normal ordering. In the boson part of 
the kinetic energy 



Tf, = ttvf / ['5p+(x) + (5p^(a;)] dx 



(5) 



the integrand has to be normal ordered. This, as above, 
is usually suppressed. 

Inhomogeneous Luttinger liquids can be realized by 
adding an external one-particle potential or by allow- 
ing the two-body interaction to depend on both spa- 



tial variables, i.e. v{x — x') in Eq. |^ is replaced by 
v{x,x'). The operator of a one-particle potential can 
only be expressed in terms of the Spa if it is suffi- 
ciently smooth in real space, i.e. the external potential 
has a vanishing 2fci? Fourier component. Similarly only 
for sufficiently smooth variations — in Sec. IIII Bl we will 
specify the meaning of "sufhciently smooth" — of v with 
{x + x')/2 the analogous steps from Eq. QJ to Eq. Q 
i.e. gi{x — x') — > gi{x,x') are allowed without changing 
the low energy physics. The standard LLL model is ob- 
tained by assuming gi{x, x') — gi^xS{x — x') 



J 



H, 



LLL,b 



vp + 1^) [5p\{x) + 5pl{x)\ + i^5p^{x)5pA^)] 



dx 




d<^^-\x) 



dx 



dx 



(6) 

(7) 



where we have defined <i>'^^'(x) = $_|_(x)±<l'_(a;), as well 
as the velocities vn{x) and vj{x). As in the homoge- 
neous model the right and left movers are coupled by the 
52 interaction and using the $'^^)(x) as the basic fields 
simplifies the solution as $'^'^^(x) and ^^'^\x') commute. 
Apart from a term which vanishes in the thermodynamic 
limit <i>^~'')(x') and d^^"'' (x) / dx obey canonical commu- 
tation relations after proper normalization 



ox 



Sl{x 



(8) 



This follows from Eqs. Q and Q. Here 6l denotes the 
L-periodic delta function. Neglecting the correction term 
yields the following Heisenberg equations of motion for 
the ^'■''^x,t) 






|*(^)(.,.) = 



-VN^x) 



d<^'-+\x,t) 
dx 



(9) 



Therefore, the ^^'^\x,t) obey generalized wave equa- 
tions, e.g. for the field related to the change of the total 
density 

^*<«(.,<,-,,(.)^.„w2i^ = o. (10) 

The spatial derivative of the first equation in Eq. @ 
yields the continuity equation for the total charge density 
6p = Sp+ + Sp^ 



9 , r .^ 9 



vj{x) 



a$(-)(a:,t) 
dx 



= 



(11) 
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FIG. 1: (Color online) Spatial dependence of velocity 'i;jv(3;) 
used for the discussion of the generalized wave equation 110^ . 
Two different realizations (full and dashed curves) of the tran- 
sition region [a, b] are shown. The shaded area shows the 
incoming density 5p{x,0). 



which implies the conservation of the total charge Q = 
J Spdx. As Eqs. (|10|l and (|ll|l are linear the expectation 
values discussed in the following obey the same equations. 



In a homogeneous system the velocities vj and wjv are 
constant and the corresponding "sound velocity" called 
the charge velocity Wc is given by Vc — ^/vjvn- It con- 
stitutes the first LL parameter characterizing the sys- 
te m. Th e second LL parameter K is defined as K = 
\/vj/vn = Vc/vjs^. The general solution of the wave 
equation for constant Vc and K reads f{x — Vct) -I- ^(x -I- 
Vct), where / and g are arbitrary functions. 



B. Computing the transmitted charge 

For general dependencies of vj(x) and vn{x) on the 
position the generahzed wave equation H1U|) cannot be 
solved analytically. In the following we discuss proper- 
ties of the solution in the thermodynamic limit when the 
velocities vj{x) and V]^{x) are constant for x < a and 
X > b > a but have arbitrary (bounded) variation in the 
interval [a, 6]. Figure ^ shows two qualitatively different 
types of behavior of V]s[{x). The dashed curve presents 
a monotonous transition between the asymptotic values 
while the full curve corresponds to an "interacting wire" 
region in [a, b] with a higher (constant) value of vn,w 
corresponding to a stronger repulsive interaction. In or- 
der to describe the charge transport through this region 
we consider an incoming density with compact support 
which at i = is completely to the left of a and moves 
towards it with constant velocity v^^l (see Fig-QJ. Before 
the disturbance reaches point a the time evolution is 

<^^ + \x,t) ^ fin{x - Vc.Lt) , Sp{x,t) ^ f!^{x - Vc.Lt) ■ 

' (12) 
The total incoming charge Qin is given by Qin = /in (a) — 
/in(— oo). This short time solution for $(+)(x,t) deter- 
mines $^~'(x,t) up to a constant. In the long time limit 
there will again be no charge in the region [a, 6] with the 
total transmitted charge Qtrans moving to the right with 
constant velocity Vcr and the total reflected charge Qrcf 
moving to the left with velocity Vq.l, where charge con- 
servation implies Qtrans + <9ref = Qin- In the following 
we show that it is possible to determine Qrct/Qin without 
the explicit solution for 6p{x, t) in this long time limit and 
without any additional assumptions on the spatial vari- 
ation of Vf^(x) and vj{x). The trick is to introduce the 
quantity 



S{t) 



1 (9$(+)(a;,i) 



dx 



-oo vjix) dt 

$(")(-cx),t)-$(")(oo,i) 



(13) 



and to show that S is time independent. This can be 
seen by discussing the expression in the second line which 
follows using Eq. (Q or if one wants to argue with the 
density only by taking the time derivative and using the 
generalized wave equation 



S{t) 



1 a2^(+)(a;,^) 



dx 



d 



-oo Vjix) 

d_ 

, dx 
vn,rSp{oo, t) ~ vn,lSp{-oo, t) ^ , 



(14) 



vn{x) — ^^+Hx) 



dx 



where we have used that the assumption about the initial 
state implies that p{±oo,t) vanishes for all finite times. 



The constant value of S follows from Eq. (|12|l 

S = - -^f[ni^-Vc,Lt)dx 

J-oo 'VJ,L 

- ^c,L ^ _ Qin 

— Vin — „ 
V.J,L Kl 



(15) 



When the incoming density hits the transition region 
[a, b] it is partially reflected and transmitted, where the 
detailed behavior is very different for the two realizations 
of Wat (x) shown in Fig.^ For sufficiently large times the 
total charge in the transition region goes to zero and 
$^+-'(a;,i) takes the form 



■I>W(a;,t) = 



const. , 

/trans (a; - Vc^rI) 



for x < a 
for a < X < b 
for X > b . 



(16) 



With the definition of 5* in Eq. (fT!^ this yields 



S 



Vc.L 
VJ.L 



Vc.,R 
VJ,R 


[ 


V; 


Qicf 


^;trans 



f^^f{x + Vc,Lt)dx 

{x — Vc,rt)dx 



Kl Kr 



(17) 



The comparison of Eqs. H15() and H17() as well as charge 
conservation leads to the result derived earlier assuming 
a stepwise variation of the LL parameters^^-^'^jii 



Qrcf 



Kl-Ke 
Kr + Kr. 



^;in 7 ^/t 



2Kf 



Kr + Kr 



Qin. (18) 



Our result shows that the properties of the transition 
region play no role at all. The ratio Qtrans/Qin is positive, 
while Qrof/Qin has no definite sign. 

From this result one can easily obtain the linear con- 
ductance G through the system. We first discuss a ho- 
mogeneous system with constant LL parameters K^ and 
Vc,L (the use of the index L becomes clear later). We 
consider a current free initial state in which the density 
is increased by 5pQ in the left half of the infinite system. 
Then half of the additional density moves to the right 
with velocity Vc.l which corresponds to the current 



i = 7:'"c,lSpo 



:Vc,L-^ Spo = G Spo 



'dp 



(19) 



in the central region extending linearly with time. Here p 
denotes the chemical potential. With dp/dp = 1/{t:vn) 
this yields for the homogeneous system 



^hom 



Kl_ 

2tt 



(20) 



We now switch to an inhomogeneous system as in Fig. ^ 
In analogy to the initial condition in the homogeneous 
case we raise the density for x < a by 5pQ. The stationary 
current is then obtained by multiplying the result for the 



homogeneous case by the fraction of transmitted charge 
through [a,b], computed in Eq. ((TH|l . For the conductance 
this leads to 



Gi 



inhom 



Kl +Kr 27r ' 



(21) 



which again is independent of the details of the transition 
region. 

For an interacting quantum wire in region [a, b] at- 
tached to noninteracting leads [Kl = Kr = 1) the con- 
ductance is l/(27r) independently of how quickly the LL 
parameters vary spatially near the contact points a and 
b. This shows that the LLL description cannot produce 
one-particle backscattering from the contact region gen- 
erated by an inhomogeneous two-particle interaction. As 
we will discuss next using a microscopic model and the 
fRG approach contacts defined by a vanishing interaction 
to the left and a positive finite interaction to the right (or 
vice versa) generically produce an effective single-particle 
backscattering. This clearly shows the limitations of the 
LLL picture. With the approximation to describe the 
two-body interaction as a quadratic form in the bose- 
fields the dependence of the conductance on the sharp- 
ness of the transition is lost. In contrast, our microscopic 
model directly allows to study the transition from smooth 
to abrupt contacts with the concomitant change of the 
conductance. 



III. A MICROSCOPIC MODEL 

As our microscopic model we consider the spinless 
tight-binding Hamiltonian with nearest-neighbor hop- 
ping r > and a spatially dependent nearest-neighbor 



interaction U, 



ij+i 



= [/. 



j+ij 



^ = - E (4+1' 

j=-oo 

+ Y, U,.j+i (n, - 1/2) (n,+i - 1/2) , (22) 
j=i 

where we used standard second-quantized notation with 
c, and c being creation and annihilation operators on 

site j, respectively, and rij — d c the local density oper- 
ator. The interaction acts only between the bonds of the 
sites 1 to A'', which define the interacting wire. We later 
take 



c/,. 



t//i(j), j = l,2,...,Af-l 



(23) 



with a function h that is different from 1 only in regions 
close to the contacts at sites 1 and A^. We here mainly 
consider the half-filled band case n = 1/2. In the inter- 
acting region the fermions have higher energy compared 
to the leads. To avoid that the interacting wire is de- 
pleted (implying a vanishing conductance) we added a 



compensating single-particle term. It can be included as 
a shift of the local density operator by 1/2. This en- 
sures that the average filling in the leads and the wire is 
ri = 1/2. It is important to have a definite filling in the 
interacting wire as the LL parameters depend on n. For 
general fillings the required shift of the density depends 
on n and the interaction.iS. 

The model Eq. H22() with nearest neighbor interaction 
U across all bonds (not only the once within 1 and A^) 
is a LL for all U and fillings, except at half filling for 
\U\ > 2t3 At n = 1/2 the LL parameter K is given 



K 



2 ( U 

— arccos 

IT V 2t 



for \U\ < 2t. To leading order in U/t this gives 



U 



K = l +0 {[U/tY 



(24) 



(25) 



At temperature T = the linear conductance of the 
system described by Eq. I|22(l can be written as2i 



G{N) = l-\t{0,N)\' 



(26) 



with the effective transmission \t{e,N)\'^ = (4r^ — [e -I- 
m]^)|Si,Jv(£, N)\'^. The one-particle Green function Q has 
to be computed in the presence of interaction and in 
contrast to the noninteracting case acquires an A^ de- 
pendence. To compute G and thus G we use a recently 
developed fRG scheme. ^^ The starting point is an exact 
hierarchy of differential flow equations for the self-energy 
S^ and higher order vertex functions, where A e (oo,0] 
denotes an infrared energy cutoff which is the flow pa- 
rameter. We here introduce A as a cutoff for the Mat- 
subara frequency. A detailed account of the method is 
given in Ref-jT^- We truncate the hierarchy by neglecting 
the fiow of the two-particle vertex only considering S^, 
which is then energy independent. The self-energy S "^ 
at the end of the fRG flow provides an approximation for 
the exact E. This approximation scheme and variants of 
it were successfully used to study a variety of transport 
problems through Id wires of correlated electrons. In 
particular, in all cases of inhomogeneous LLs studied the 
exponents of power-law scaling were reproduced correctly 
to leading order in [/J£i24j25j26 

On the present level of approximation the fRG flow 
equation for the self-energy reads 



aA^^'-i" 2^ ^.^ 



e b'2 2'(«w)ri/,2';l,2 , 



L)=±A 2,2' 



(27) 
where the lower indices 1, 2, etc. label single-particle 
states, ri/^2';i,2 is the bare antisymmetrized two-particle 
interaction and 



G\i^)^{Q^\iuj)-T,'^\ ' , 



(28) 



with the noninteracting propagator Qq. For the ini- 
tial condition of the self-energy flow see below. As our 
single-particle basis we later use Wannier states [as in the 
Hamiltonian 122|l ] as well as momentum states. 



Numerical results 



In the real space Wannier basis the self-energy matrix 
is tridiagonal and the set of coupled differential equations 
(jTTjl reads (with j € [l,iV]) 



dA 



2tt 
dA ^'^^^ 2tt 



E ^i^,±iMe'"°"-(30) 



w=±A 



To derive these equations we have used that 

with C/jij2 = C/jj jj_|_i(i5j^_J2_i -I- 6j-^j2+i)- The initial 
condition of S^ at A = oo isifliS^ 



^3d - 



(Uj- 



ij 



U. 



jj+i 



)/2 



"ij±i 



= 



(32) 



Even for very large TV (earlier results were obtained for 
up to A^ = 10'^)i2i^ the system Eqs. ^ and (EOl can 
easily be integrated numerically starting at a large but 
finite Aq. Due to the slow decay of the right-hand-side 
(rhs) of the flow equation the integration from oo to Aq 
gives a nonvanishing contribution on the diagonal of the 
self-energy matrix that does not vanish even for Aq — > oo 
and that exactly cancels the term in Eq. H32(l such that 






"JJ±1 



0, 
0. 



(33) 



For half filling the Hamiltonian (|22() is particle-hole 
symmetric. For N odd and symmetric h{j), that is 
h{j) = h{N — 1—j), on which we focus here, it is further- 
more invariant under inversion at site j ~ (A^ + f )/2. To- 
gether these symmetries lead to a vanishing conductance 
as was discussed earlier i2L^ We thus only consider even 
A^. The two extreme cases of (a) an abrupt turning on 
and off of the interaction h{j) = 1 for j = 1,2, ... ,N —I, 
and zero otherwise, and (b) a very smooth variation of 
h{j) were considered earlier. In case (a) l/(27r) — G in- 
creases for fixed A^ and increasing^^-^* U as well as for 
fixed U and increasing A^.^^ Using our fRG scheme it was 
shown that for asymptotically large A^, G vanishes as 



G ^ {6n/t) 



2{1/K-1) 



(34) 



with the energy scale 6n = nvp/N. The scale on which 
the asymptotic low-energy scaling sets in strongly de- 
pends on U and even for up to A^ = 10^ it was only 
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FIG. 2: (Color online) Upper panel: The conductance 1 — 
2ttG as a function of Sn/t for U/t — 0.2, contacts of ?n = 
7 lattice sites, and three different contact functions hi{j) of 
increasing smoothness. Lower panel: The effective exponent 
(logarithmic derivative) of the data. The thin solid line shows 
the exponent 2{K - 1)/{K + 1) with K{U/t = 0.2) = 0.9401. 



reached for fairly large U, e.g. U = 1.5t. As discussed 
in the introduction scaling of G with the same exponent 
as in Eq. (|Mjl is found for a single impurity in an oth- 
erwise perfect chain. This suggests that the low-energy 
physics of a correlated quantum wire with two contact in- 
homogeneities due to the two-particle interaction is simi- 
lar to that of a wire with a single impurity (single-particle 
term). Here we will further investigate this relation. 

Already at this stage it is important to note that the 
similarity is not complete. As shown elsewhere^^ G{T) 
for two contacts, fixed U, fixed large A^, and decreasing 
temperature T vanishes as (for T '^ Sn) 



G - (T/t)1/^~i 



(35) 



that is with an exponent that is only one-half of the one 
found in the Sn scaling Eq. (|34|l . This is in clear contrast 
to the temperature dependence of G for a single impurity 
in an otherwise perfect chain. In this case the exponent 
2(1/Ar— 1) is found in the scaling with any infrared energy 
scale, e.g. T, 6n, and A. 

For a given A and U it is always possible to find a 
sufficiently smooth function h{j) such that f/(27r) — G 
is smaller then an arbitrarily small upper bound [case 
(h)]i^i^ In that case the contacts are perfect and the 
effect of impurities in the interacting wire connected to 
Fermi liquid leads can be studied. Further down we will 
give a quantitative definition of the meaning of "suffi- 
ciently smooth" . 

Here we study the dependence of the conductance on 
A^, U, and the function h{j) by numerically solving the 
flow equations (|29|l and (|30|l with the initial condition 
Eq. I|33() . In Fig. [5] we show 1 — 2ttG as a function 
of 6n/t ^ 1/N on a log- log scale (upper panel) for 
U/t = 0.2 and three different h{j) of increasing smooth- 
ness. In our numerical calculations for simplicity we focus 



on contact functions h{j) that are symmetric around the 
middle of the interacting wire (equal contacts) and there- 
fore only give their definition for the first half of the wire 

(withje[i,iv-i]) 



h,{j) = 1, j = l,...,iV/2 



h2ij) = 



_ / j/m, j = 1, 



J = m, . 



, m — 1 

.,N/2 



and 



hsU) = 



2(j7m)2, j = l,...,(m-l)/2 

l-2{j/m~l)'^ , j = (to + 1)/2,...,to 
1 , j=m,...,N/2 



where m is odd and measures the length of the contact 
region. In Fig. |21 we chose m = 7. The lower panel of 
Fig. [21 shows the logarithmic derivative of the data. A 
plateau in this figure corresponds to power-law scaling of 
the original data with an exponent given by the plateau 
value. For fixed Sn, ^ — iiiG decreases quickly with in- 
creasing smoothness of h{j). For sufficiently small 5n /t 
all curves follow power-law scaling. The smoother the 
contact the smaller 5m /t has to be before the scaling sets 
in. The numerical exponent is close to 2{K — 1)/{K -1-1) 
(shown as the thin solid line), with K from Eq. (|24|l . 
The latter is the exponent found for a system with per- 
fect contacts and a weak single-particle impurity placed 
in the contact region. We also studied other values for 
the length of contact m and found similar results. The 
deviation 1 — 27rG from the unitary limit decreases if m 
is increased while all other parameters are fixed. 

Within our approximation scheme we map the many- 
body problem on an effective single-particle problem 
with the energy independent S as an impurity 

potential. ^""-^ During the fRG flow the interplay of the 
two-particle inhomogeneity and the bulk interaction gen- 
erates an oscillating self-energy with an amplitude that 
decays slowly away from the two contacts. This is in close 
analogy to the case of a single impurity in an otherwise 
perfect chain»ifiiSiSl Scattering off this effective potential 
leads to the reduced conductance. 

The U dependence of 1 — 2tiG and the effective expo- 
nent is shown in Fig.|31for /i2(i). At fixed 5n the devia- 
tion of the conductance from the unitary limit increases 
with U . For all U /t shown we find power-law behav- 
ior with exponents that are close to 2{K — l)/(^ +1) 
(again shown as the thin solid lines). The larger U/t 
the larger is the deviation of the numerical exponent and 
2{K — 1)/{K +1). This is not surprising as our approxi- 
mate scheme can only be expected to capture exponents 
to leading order in U /t}^ 

For later reference we note that power-law behavior is 
not only found in the Sm (that is 1/A^) dependence of 
1 — 2nG but also in the dependence on the fRG flow pa- 
rameter A at a large fixed N . In this case G is computed 
using Eq. H2ti|) with the Green function at scale A. This 
is shown in Fig. 0] where we compare the dependence on 
Sn/t and A/r for U = 0.2 and ft.2(j)- Scaling with A/r 
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holds roughly down to the scale 5n/t = itvf/{tN) set 
by the length of the interacting wire. Beyond this scale 
G saturates. 

We now analytically confirm the power-law scaling of 
l/(27r) — G for weak effective inhomogeneities and com- 
pute the exponent to leading order in the bulk interac- 
tion. 



B. Analytical results 

Numerically we found that l/(27r) — G shows power- 
law scaling in both effective low-energy cutoffs 5^ and A, 
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as long as the scaling variable considered is sufficiently 
larger than the other fixed energy scale. We note that the 
same holds for the temperature T as the variablei^ In all 
cases the same exponent is found. Analytically the easi- 
est access to scaling is by considering A^ ^ oo, T = 0, and 
taking A as the variable. We will here follow this route. 
For the case of a single impurity with bare backscattering 
amplitude Vback in an infinite LL (no Fermi liquid leads) 
a similar calculation was performed within the fRG^^ to 
show that for small A 

SA..,-fe./r^ (Fback/r) (A/t)-^/(--) . 

Within the Born approximation this leads to 

f -G^|yback/r|'(A/r)-2^A-), 

where — 2C//(7rT) is the leading order approximation of 
the weak impurity exponent 2{K— 1) [see Eq. (|25|l ]. This 
scaling holds as long as the rhs stays small. For analytical 
calculations it is advantageous to switch from the real 
space basis to the momentum states. 

In the momentum state basis the flow equation 1)27(1 is 
given by 



(36) 



(37) 









j=±A 



J „iwO+ r'A 



dgdg'e^"*^ g^M^)Vk'.,^';K^, 



(38) 



J{i-q') _ pi(l~k') 



with 

''k',q';k,q — „ 

xU{k + q-k' -q') 
and the Fourier transform of the interaction 



^i(k-k') _ ^i(k-q') 



(39) 



U{k) 



U h{k) . 



Jjk 



(40) 



Here h{k) denotes the Fourier transform of the function 
h(j) that contains the shape of the turning on and off 
of the interaction. In the present section we do not as- 
sume special symmetry properties of h{j). Thus the two 
contacts might be different, as it is generically the case 
in experiments. At A = Aq all matrix elements of S^" 
are zero. We now expand the rhs of Eq. (|38|l to first 
order in E^. This expansion is justified as long as T,^ 
remains small (it is certainly small at the beginning of 
the fRG flow). The flow equation then becomes an inho- 
mogeneous, linear differential equation 

^-' 4}1{A) + 4%{a,j:^) , (41) 



yA 



with 



^i'^A) = - 



-E 

27r ^^. 



^M.-.')(e- 



^ — ik 



w = ±A 



1- 



dq 



k' ,q-^k,q 



lUJ - 



i. 



A 



VA2 + 4t2 



,(42) 



where we used ^k — — 2TC0s(fc) — ^ and Eq. l(^ . The 
factor e*'^° was dropped as the integration starts at 
Ao < oo. Later we will primarily be interested in the 
kp,—kp matrix element (backscattering) of S. For these 
momenta the rhs simplifies to 



^E!-..(A) 



U r 



h{-2kF)e 



— ikf 



A 



VA2 + 4r2 



(.43) 



The term on the rhs of the linearized flow equation that 
contains S'^ is 



•^S(A,S^) = "^ E r dqdq'rk,,q,.,k^q 



(SoS^So),,,, 
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(44) 



Before further evaluating this term we take Af — > oo and 
compute T^i ^q-^k,q in this limit. As we will see in this case 
the information about the smoothness of h{j) is lost on 
the rhs of Eq. (I44|) . It nevertheless enters the solution of 

the differential equation via the inhomogeneity ^k\ in 
Eq. H41() . As we are only interested in the leading diver- 
gent behavior of the linear part on the rhs of Eq. H41|l 
this procedure is justified. To be specific we first assume 
that the interaction is turned on and off abruptly, that is 
h{j) = 1 for j = 1, 2, . . . , A^ — 1 and zero otherwise. To 
obtain the Fourier transform of the two-particle interac- 
tion we have to perform 



Uik) 



Uh{k) 



U 

2^ 



JV-l 

E' 



^ijk 



U e"" - e 



\Nk 






27r 1 



(45) 



With this result the two-particle vertex simplifies to 
Ffc 



k',q;k,q "^ — [cOs(fc - fc') - CQs{q - k')] 



X(527r(fc + q — k' — q) 



(46) 



Up to a factor 1/2 the vertex is then equivalent to the one 
obtained for a homogeneous nearest neighbor interaction. 

With this vertex J^^?^^^.(A, S^) Eq. ^ reads 

:FiJ{(A,j:^)^-^J2 r dq[cos{k-k') 



t^=±A- 



cos(g — k')] 



-'q+k — k' ,q 



{iu - ^q+k-k'){iUJ - Cr?) 



(47) 



The differential equation (|41|l can be solved by the vari- 
ation of constant method. To apply this we first deter- 
mine the solution of the homogeneous equation. The 
scaling can be extracted if only the leading singular con- 
tribution (for A ^ 0) ofTJ^^j^iA, E'^) Eq. gZJ at /c' = kp 



and k = —kp is kept. Following the same steps as in 
Ref . M we find 



^t- 



U 1 



^ s 



^-1 horn 



-fcFJhom- 2tTT a 

with the solution 



(48) 



(49) 



where cq is a dimensionless constant. The singular part 
of the solution of the inhomogeneous linear differential 
equation is then given by 



7r « c(A) (A/r)-^/(2-) , 



(50) 



with 



c(A) = / dA' 






•J 1 



(1) 



..(A'). (51) 



horn 



For small A the function c(A) is non-singular and given 
by 



lim c(A) = h(-2kF)e 



-ikp 



A— ■ ■ T 

with c being a constant of order 1 . This gives 



(52) 



^kp.-kp/''' ' 



Ur 



h{-2kF)e-''"' c (A/r)-^/(2^^) (53) 



inhomogeneity this is replaced by the square of the bulk 
interaction and the square of the 2kF Fourier component 
(backscattering) of the function h{j) with which the in- 
teraction is turned on and off. The smaller this com- 
ponent the smaller is the perturbation due to the two 
contacts. Therefore the smoothness of the turning on 
and off is directly measured by h{—2kp). The presence 
of the factor \U/t\'^ explains why in the numerical study 
the weak inhomogeneity exponent for larger U/t is only 
observable for fairly smooth contacts, that is /i(j)'s with 
small 2kp component. For larger U/t and fairly abrupt 
contacts the rhs of Eq. H54|l becomes too large already 
on intermediate energy scales and no energy window for 
scaling is left. For A <C Ac defined in Eq. H55|) the inho- 
mogeneity is effectively large and Eq. H34|) holds. 

We here considered the half-filled band case, but also 
other fillings can be studied following the same stepsi^^ 

The analysis Eqs. I|38|l to (|54|l also holds if in addi- 
tion weak single-particle impurities are placed close to 
the contacts, with the only difference that E^° _j, now 
has a non- vanishing initial condition set by the backscat- 
tering of the bare impurity. The prefactor in Eq. H54() 
is determined by either the square of the single impurity 
backscattering amplitude or \U/t\ \h{—2kp)\ depend- 
ing on the relative size. 



IV. SUMMARY 



and with the Born approximation the final result 



1 

2^ 



G-|C//r|' hi-2kp) (A/t)-^/("") . (54) 



To obtain the singular part of the solution of the homo- 
geneous differential equation S^ _^, we assumed 

L ^' -'^ J horn 

that the interaction is turned on and off abruptly. One 
can show that the parts of h{j) with a smooth variation 
of finite length do not contribute to the singular part of 
Eq. 144(1 . Thus the same singular part is found indepen- 
dently of how the interaction is varied and Eq. H54|) is 
valid for general h{j). 

The power-law scaling Eq. (|54|l holds for Ac ^ A ^ r 
with a scale Ac set by 

\U/t\^ \hi-2kp)\^ (Ac/t)-^/(-^) ~ 1/(4^) . (55) 

To leading order in U/t the exponent [//(ttt) agrees with 
the exponent 2{K-l)/{K+l) [see Eq. ^] found for a 
single weak impurity placed close to a perfect (that is ar- 
bitrarily smooth) contact i2iiS We thus have analytically 
shown, that with respect to the scaling exponent weak 
single-particle and weak two-particle inhomogeneities are 
indeed equivalent. The leading order exponent is further- 
more consistent with the numerical results of the last 
section. 

For the single impurity case the prefactor of the power 
law is given by |Vback/'''| • In case of the two-particle 



In the present paper we have investigated the role of 
contacts, defined by an inhomogeneous two-particle in- 
teraction, on the linear conductance through an inter- 
acting Id quantum wire. The wire and contacts were 
assumed to be free of any bare single-particle impurities. 
We first showed that within the LLL picture the contacts 
are always perfect, that is the conductance is l/(27r) in- 
dependent of the strength of the interaction, the length of 
the wire, and the spatial variation of the LL parameters. 
Earlier only step-wise changes of the LL parameters were 
considered. We then studied the problem within a micro- 
scopic lattice model. Similar to the case of a single impu- 
rity the interplay of the two-particle inhomogeneity and 
the bulk interaction generates an oscillating self-energy 
with an amplitude that decays slowly away from the con- 
tacts. Scattering off this effective potential leads to the 
discussed effects. We showed that within the microscopic 
modeling l/(27r) — G increases with increasing bulk in- 
teraction [/, increasing wire length N, and decreasing 
smoothness. The measure for smoothness is given by 
h(—2kp), that is the backscattering Fourier component 
of the spatial variation h(j) of the two-particle interac- 
tion. As long as the inhomogeneity stays effectively small 
l/(27r) — G shows power-law scaling with an exponent 
that is consistent with 2{K — 1)/{K -1-1), the scaling ex- 
ponent known for the case of a wire with a single impurity 
placed close to one of the two smooth contacts. 

In experiments on quantum wires the leads are elec- 
tronically two- or three-dimensional. Depending on the 
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systems studied the contacts are either regions in which 
the system gradually crosses over from higher dimensions 
to quasi Id or the contact regions extend over a finite part 
of the wire (see the experiments on carbon nanotubes)^ 
This shows that our simplified description — Id leads, spa- 
tially dependent interactions close to end contacts, no ex- 
plicit single-particle scattering at the contacts — provides 
only an additional step towards a detailed understanding 
of the role of contacts in transport through interacting 
quantum wires, such as carbon nanotubcs and cleaved 



edge overgrowth samples 
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